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FORMULAS  AND  CHARTS  FOR  THE  SUIERSONIC  LIFT  AND  DRAG  OF 
FLAT  SWEPT^CK  WINGS  WITH  INTEEACTING 
LEADING  AND  TRAILING  EDGES 
By  Doris  CoLen 


SUMMARY 


The  prohlem  is  considered  of  a  wing  with  rectilinear  plan  form  swept 
so  that  hoth  leading  and  trailing  edges  lie  within  their  respective  Ifeich 
cones;  moreover,  the  Mach  lines  from  the  traillng-edge  apex  intersect  the 
leading  edge.  Formulas  and  design  charts  are  presented  for  the  lift  in 
such  a  case,  based  on  approximate  formulas  for  the  lift  distribution 
developed  in  NACA  TN  1991^  19^9*  Ihe  charts  cover  a  practical  range  of 
aspect  ratios  and  plan  forms  of  moderate  taper,  with  tips  parallel  to  the 
stream.  The  leading-edge  thnist  and  drag  due  to  lift  are  also  readily 
calculated  from  the  material  presented.  NumericaQ.  results  ant!  an  applica¬ 
tion  of  the  charts  are  included. 


INTRODUCTION 


Range  of  Applicability  of  Formulas 


Problems  in  linearized  supersonic  wing  theoiy  are  characterized 
primarily  by  the  orientation  relative  to  the  wing  boundaries  of  the  Mach 
lines  arising  from  the  various  points  of  disturbance  defining  those 
boundaries.  Therefore,  even  when  the  plan  form  of  a  wing  is  specified, 
a  series  of  problems  arises  if  the  wing  is  to  fly  throuf^  any  consider¬ 
able  range  of  Ifech  numbers. 

Paradoxically,  the  problem  of  the  flow  at  higher  Mach  nimibers  has 
proved  more  readily  amenable  to  solution  than  that  at  lower  speeds, 
because  at  the  hi^er  speeds  the  zones  of  influence  (Mach  cones)  are 
narrow  and  interference  problems  are  fewer.  In  the  case  of  the  conven¬ 
tional  swept-bach  wing  at  an  angle  of  attack,  four  Mach  number  ranges  may 
be  said  to  have  been  investigated  to  a  sufficient  extent  that  the  lift 
and  drag  of  any  specified  wing  is  readily  obtainable.  The  aerodynamic 
characteristics  were  first  computed  (reference  l)  for  that  range  of  Mach 
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numbers  in  which  only  the  center  sections  of  the  wing  are  included  in  the 
Mach  cone  from  the  leading-edge  apex;  formulas  have  also  been  given 
(reference  2)  for  the  next  lower  range,  where  the  leading-edp  Mach  line 
Lts  the  tip;  formulas  (reference  2)  and  charts  (reference  3) 
given  for  a  third  range,  where  the  leading-^dge  Mach  line  stands  ahead  of 
the  wing,  but  the  Mach  line  from  the  tr ailing-edge  apex  is  behind  the  ^ 
wing;  Sd  formulas  have  been  given  (reference  2)  for  the  fourth  case,  in 
which  the  trailing-edge  Mach  cone  includes  part  of  the  tip,  but  not  the 
leading  edge. 

The  last  of  these  cases  is  illustrated  in  figure  l(a).  In  > 

the  same  wing  is  shown  at  a  smaller  Mach  number  such  ttet  the  trailing 
edge  Mach  line  now  intersects  the  leading  edge.  In  reference  , 
in^such  a  case  was  investigated  and  formulas  for  the  load  distribution  ^ 
were  derived.  In  the  present  report,  the  load  distribution  will  be  inte 
grated  to  obtain' formulas  for  the  total  lift.  Formulas  for  the  lea  ing 
edge  thrust  and  the  drag  due  to  lift  will  also  be  presented. 

The  Mach  number  range  to  which  the  formulas  and  charts  apply  is  indi¬ 
cated  as  a  function  of  taper  ratio  X  and  aspect  ratio  A  for  severa 
Sles  of  sweepback  A  in  figure  2.  The  upper-limit  curves  are  determined 
by  the  condition  that  the  trailing-edge  Mach  line  intersects  the  wing  P 

at  the  leading  edge.  Above  these  curves,  the  formulas  of  reference  2  y 

be  used.  A  lower  Mach  nuniber  range,  merging  into  the  transonic  ®uc 

that  the  flow  field  of  one  wing  tip  extends  laterally  to  include  a  part  o 
Sr oSer  tip,  is  defined  by  the  lower-limit  curves.  The  flow^pattern  in 
such  a  case  is  too  complicated  to  be  treated  by  tte  method  of 
report.  Mach  numbers  low  enough  to  permit  wing-tip  interference 
within  the  scope  of  the  approximate  theory  of  reference  5- 

The  problem  of  wing— tip  interaction  also,  of  course,  sets  a  lower 
limit  m  lit  range  of  Spec?  ratios  covered.  On  the  other  h^,  the  me^od 
of  reference  4  and  the  present  paper  may  in  theory  be  applied  to  the  com¬ 
plete  range  of  aspect  ratios  to  the  right  of  the  boundaries  shown  in 
ficure  2.  However,  the  necessary  fommlas  have  been  developed  only  up 
a  certain  point  beyond  which  the  services  of  high-speed  coi^uting  machin 
ery  would  be  required  for  their  evaluation.  Both  mathematical  and  physical 
considerations  make  the  application  of  the  method  to  very  highly  tapered 
plan  forms  inadvisable.  For  cases  of  moderate  taper  ^d  a  practical  r^e 
of  aspect  ratios,  charts  have  been  computed  to  facilitate  the  calculation 
of  both  lift  and  drag  due  to  lift.  The  formulas  and 

developed  for  wings  with  tips  parallel  to  the  stream  only,  but  the  general 
procedure  is  applicable  to  wings  with  raked  tips  also. 

The  s3niibols  used  are  listed  in  the  appendix. 


NACA  OT  2093 


3 


General  Outline  of  the  Procediire 


It  ^  interaction  takes  place 

fSr®''  *  leading  and  trailing  edges,  the  ving  plan 

vr  co^rise  tvo  principal  regions,  separated  (fig.  1(h))  by 

Ifec  line  arising  at  the  point  of  intersection  of  the  trailing-edge 
Mach  l^e  aM  the  leading  edge.  Ahead  of  this  line  (region  l)  tS  flL 

(reSon  described  in  terms  of  conical  fields.  Behind  this  line 

iregion  ip,  the  flow  is  more  nearly  two-dimensional,  and  the  loading  can 
in  the  well-known  subsonic  flat-plate  formula,  corrected 

in  ma^itude  to  give  the  proper  value  of  the  pressure  at  the  leading  edge 
Of  each  section  The  formulas  for  correcting  the  two-dimensional  ?Ld£ 

SafirL'"''  “  i^iaitely  long  wing,  of  course,  the 

g  pproaches  the  subsonic  flat— plate  loading,  corrected  for  the  angle 
of  yaw  according  to  simple  sweep  theory.  ^ 

loading  in  region  I  the  conical- 

flow  solutions  of  reference  2  and  in  region  II  the  formulas  given  in 
reference  4.  From  the  integrated  lift  and  the  leading-edge  thrust,  the 
^ag  due  to  lift  can  be  calculated.  Hhe  leading-edge  thrust  in  region  I  is 
the  sa^  as  that  for  the  tri^gular  wing  with  the  same  leading-^dgf  sweep' 
^d  has  been  given  by  Jones  (reference  6),  Hayes  (reference  7),  Eobinson 
(cited  in  reference  8),  and  others.  The  disturbance  arising  at  the 
trailing  edge  affects  the  leading-edge  suction  in  region  II.  Formulas  are 

calculating  the  resulting  thrust,  from  which  the  drag  due  to 
lift  follows . 


FORMULAS  FOR  LIFT 


Region  I 


.  .  flovs  method.-  The  term  "conical-flows  method”  is  used  to 

T®  briefly  a  method  of  superposition  of  conical  flow  fields ^  devel- 
Lagerstrom  (reference  9)  as  an  extension  of  some  earlier  work  by 

is  distribution  on  a  swept-back  wing 

iL  ^  superposing  on  the  conical  loading  of  an  infinite  triang^- 

ll  c^el  flow  fields  of  such  magnitude  and  orientation  as 

0  c^cel  loading  on  the  portions  of  the  triangular  plate  lying  outside  the 
boundaries  of  the  specified  swept-back  wing.  ?his  process,  as  apSiefaJ 

reference's!  ^  swept-back  wing,  is  more  fully  described  in 

^A  conical  field  is  one  in  which  the  velocities  are  constant  alone  ravs - 

emanating  from  a  point  (the  apex  of  the  cone). 


k 


NACA  TM  2093 


In  figure  l(t),  region  I  has  been  further  divided  into  regions  la 
and  Ib  hy  the  Mach  line  from  the  trailing-edge  apex.  Eegion  Ib  is  t 

2  Leu  the  loedlBg  le  .^.difled  by  the  addition  of  tields^ed 
to  cancel  the  load  behind  the  trailing  edge.  In  region  la,  ^ 
siii5)ly  that  on  a  triangular  wing  having  the  same  apex  angle  as  the  swept 

back  wing.  .  . 

The  conical  flow  fields  were  originally  derived  «st  ccnvenl^t^  M 
velocity,  rather  than  pressure,  fields.  The  loading  Is  obtained  from  the 
velocity  by  the  linearized  relation 


Ap  _  Au 
q  "  V 

where  u  is  the  streafflwise  component  of  the 

upper  surface,  V  is  the  stream  velocity,  and  Ap/q  is  the  coeffici 
of  lifting  pressure. 

T.oad  distribution  in  region  I.-  In  terms  of  the  conical  variable 
a  =  — ilie  loadi^  over  tne  uriangular  wing  is  given  by 


mu. 


^A  ,-^.2 


(1) 


where 


Ur.  = 


mV^a 


(2) 


pE( 


is  the  strearwise  coii5)onent  of  the  perturbation 

center  line  (a=0)  and  u^  is  the  same  velocity  component  elsewhere  on 
the  wing. 

Cancellation  of  the  load  in  the  wake,  as  described  in  reference  2, 
is  accon^jlished  by  superposing  first  a  symmetrical  field  teving  ®  ' 

load  reSi  coincident  with  the  wake  to  cancel  the  perturbation  velocity 
u«  throughout  the  wake,  and  then  an  infinite  number 
Saded  oblique  elements  to  c^cel  the 

constant-load  region  (shaded)  of  a  single  oblique  flow  field. 

Cancellation  of  Uo  induces  on  the  wing  a  press^e  distribution 
(equation  (51)  of  reference  2,  corrected)  proportional  to 


(Au)_  =  - 


u. 


K(  V 


Y  (  >/  cp) 


(3) 
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where  K(  is  the  complete  elliptic  integral  of  the  first  kind 

of  modulus  V 1— mt  ,  and  F  is  the  incon^jlete  integral  with  argument 


(p  =  sin 


Each  of  the  oblique  elements  Induces  on  the  wing  the  perturbation 
velocity  (equation  (50)  of  reference  2) 

(Au),  =  ^  cos-1 

”  ( 1-114.  )(tg^-a) 

where  Uq  is  the  differential  of  u^  obtained  by  differentiating  with 
respect  to  or 


“a  - 

Uncorrected  lift  in  region  I.-  In  order  to  obtain  the  total  lift  in 
region  1,2  we  first  integrate  the  uncorrected  loading  of  the  triangular 
wing  over  the  entire  area,  obtaining 


L_  If-m  Uq  2  f  1+m  1 

-id.  - - ^  Cq  V  -  .  ^  pQg 

qa  p  Va  L(l-ffl)®  L  V  i-m^ 


2  2 
m^-€L2 


m( l+a2 )  1+^2 


^“6-2 


X- 


2  ^2 
mt  -HOI 


m  — m^as  ,  m 

- :  —  COS  - 

m(  a^)  ^ 


where  as  is  the  value  of  a  corresponding  to  the  ray  through  the  point 
of  intersection  x^^ys  of  the  trailing  edge  and  the  Mach  line  separating 
regions  I  and  II.  (See  figs.  1  and  3»)*  In  terms  of  the  wing  plan— form 
parameters , 


a.2  = 


2mii^. 

l+m+jo^—cim^ 


formulas  which  follow  are  for  the  complete  wing;  that  is,  the  left^iand 
side  is  included. 
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When  iat;=m  (untapered  wing),  the  second  part  of  equation  (j)  becomes 
indeterminate.  In  this  case. 


hi 

^  2. 

f  l+m 

1 

2 

—1  ^  ’^Q'2  ' 

/  2  2  “1 
V  ni  -€L2 

■■■Vi  zr 

qa 

P  Va  ° 

Ul-m)3 

^  V 1— 

cub  . 

ni(  1+3,2  / 

l+a2 

+ 


1 

3 


m 


m-aa 


(9) 


t 

Wake  corrections.—  The  wake  corrections  are  to  be  inte^ated  over 
region  lb.  Integration  of  the  symmetrical  wake  correction  (eqmtion  (3)) 
yields 


(ADq  ^  -l6m^3:i^  ^  L 
qa  p (1+m^)  Va  \ 


E  (  'i/L— 

K(  J  f 


where 


(10) 


(11) 


is  the  X  coordinate  of  the  intersection  of  the  trailing-edge  Mach  line 
with  the  leading  edge.  (See  fig.  3*) 

For  each  oblique  element,  the  reduction  in  lift  is  given  by 


i  Mi 

qa  da 


da 


in  which 


gua  l+^t 

pVa  1+a 


(mt-a)(x2-xa) 


(l+mt)Cl-a) 

2(mt-a) 


(12) 


(13) 


is  the  X  coordinate  of  the  intersection  of  the  Mach  line  from  Xi,yi 
with  the  trailing  edge  (see  fig.  3)  and 


Xa 


is  the  X  coordinate  of  the  apex  of  +be  element  a. 


(14) 
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Lift  in  region  I.-  The  total  lift  in  region  I  is  then  given  by 


2  dAL 

+  —  /  - da 

qa  Jq  da 


The  quantity 


^ — -f — ^  is  plotted  against  m-h  in  figiire  4  for 


several  values  of  the  ratio  This  latter  parameter  is  the  ratio 

of  the  tangents  of  the  semiapex  angles  of  the  leading  and  trailing  edges 
and  is  constant  for  any  one  wing  through  the  Mach  number  range . 


Region  II 


Tip-effect  prohlems,—  In  figure  1(1))^  region  II  is  shown  divided 
into  two  sections  hy  the  Mach  line  from  the  leading  edge  of  the  tip. 

Such  a  division  is  actually  an  oversimplification  of  the  problem.  While 
region  Ila  can  be  treated,  by  the  method  of  reference  4,  in  such  a  way 
as  automatically  to  satisfy  the  Kutta  condition  at  the  trailing  edge,  the 
tip  effect,  which  modifies  the  lift  in  region  Ilb,  has  not  been  determined 
so  as  to  take  this  condition  into  account.  In  the  cancellation  method 
used  to  determine  the  tip  effect,  a  so-called  "primary  tip  correction"  is 
first  obtained  by  superposing  conical  flows  to  reduce  the  lift  to  zero 
along  and  outboard  of  the  side  edge.  A 
further  succession  of  steps  (reference  2) 
is  required  to  cancel  lift  introduced  in 
the  wake  behind  region  lib,  and  thereafter 
outboard  of  the  tip  within  the  Mach  cone 
from  x*,y*  (see  sketch)  and  so  on.  The 
corrections  obtained  by  this  procedure 
alternate  in  sign  and  become  successively 
smaller  in  magnitude,  while  increasing  in 
mathematical  complexity.  Since,  moreover, 
it  is  known  from  experiments  that  the 
assumed  flow  in  the  tip  regions  is  at  var¬ 
iance  with  the  physical  flow,  it  would  be 
illogical  to  attempt  any  precise  evalua¬ 
tion  of  these  corrections*  A  simple  for¬ 
mula  will  therefore  be  given,  following 
the  derivation  of  the  primary  tip  correc¬ 
tion,  for  obtaining  a  fair  estimate  of  the 
"secondary  correction"j  further  correc¬ 
tions  will  be  neglected. 

Loading  in  region  ll.—  Except  for  tip 
losses,  the  lift  in  region  II  may  be 
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calculates  ty  iutegratlng  out  to  the  tip  the  approxluate  load  distribution 
given  in  reference  4.  This  loading  is 


^  [Py-^(^-Co)lCo 

[  ( x-c  Q )  ]  ( mx-py ) 


which  is  merely  the  subsonic  two-dimensional  flat-plate  loading,  applied 
on  sections  taken  normal  to  the  stream,  and  with  a  correction  factor 
inserted  to  bring  the  loading  at  the  leading  edge  into  agreeme^  with  that 
derived  for  the  swept-^aack  wing  by  the  conical-flows  method.  The  theoreti 
cal  loading  at  a  subsonic  leading  edge  is,  of  course,  infinite  in  w  case, 
with  the  infinity  entering  as  the  reciprocal  square  root  of  the  dis-^ce 
to  the  leading  edge.  The  coefficient  or  is  a  measure  of  the  strength  of 
the  leading-edge  singularity  in  the  loading  on  the  swept-back  wing. 

It  is  most  convenient  to  describe  the  leading-edge  singularity 
terms  of  the  coefficient  in  the  perturbation  velocity  u  of  ^ 

where  the  quantity  (mx-py)  is  p  times  the  s£a^ise  distance  to  the 
leading  edge.  This  coefficient  is  then  u-/m-py  and  the  stren^h  of 
the  singularity  in  u  is  the  limit  of  its  value  as  py  approaches  m 

(the  leading  edge).  The  coefficient  a  is  a  nondimens ional  form  of  the 

same  quantity: 


iax-t>y 

Cq  V 


and  is  a  function  of  x. 


py=inx 


In  line  vith  the  general  procedure  of  the  conical— flows  method^  the 
coefficient  Ca  of  (rax-Py)“^/^  corresponding  to  the  triangular-wing 
velocity  distribution  is  found  first,  and  is  then  corrected  to  take 
account  of  the  effects  of  canceling  the  loading  in  the  wake,  by  the  sym¬ 
metrical  canceling  element  and  by  the  oblique  elements .  The  resulting 
expression  for  cr  is 


c(x)  =  - ^  +  (ZsC)o  + 

Vex vCq  L 


®-o 


where 


=  ^ 


(19) 
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(see  sketcli)  and  Turing  about  further 
modifications  of  the  flow  on  the 
outboard  portions  of  the  wing.  The 
only  effect  of  this  modification  on 
the  formulas  presented  herein  is  the 
introduction  of  additional  tenns 
into  the  expression  for  a  for^ 
values  of  x  beyond  X3,  and  in 
still  more  extreme  cases,  beyond 
X5,  and  so  on.  Evaluation  of  these 
terms  by  the  conical-flows  method 
involves  multiple  integration  and  is 
impractical  except  with  the  aid  of 
high-speed  computing  machinery. 

Each  successive  coirection  (to  c) 
is  initially  zero  and  enters  with 
zero  slope  at  X3,  zero  slope  and 
curvature  at  x  5  and  s  o  on,  so  that 
the  three-term  expression  for  cr 
given  by  equation  (18)  may  be  used 
with  satisfactory  accuracy  for  some 
distance  beyond  the  last  value  of 
X  (xg)  for  which  it  is  strictly 
valid.  In  practice,  the  third  term 
in  equation  (I8)  may  also  be  neg¬ 
lected  for  values  of  x  only  slightly 
greater  than  xi. 


Charts  giving  ^  as  a  function  of  for  the  values  of 

the  tapei^ratio  parameter  —  covered  in  figure  4  are  presented  in 

figure  5  as  an  aid  to  conqputing.  The  curves  vere  computed  using  equa¬ 
tion  (18)  and  are  therefore  exact  only  up  to  x=X3  (shown  by  a  vertical 
mark  in  each  curve).  The  points  x=x  are  also  indicated  (by  X*s)  as  a 
more  practical  limit  to  which  use  of  the  curves  may  be  extended.  (These 
points  are  off  the  scale  for  mt  =  0,8  and  0.9  figure  5(a)).  When  the 
wings  are  untapered  (m/m^-  =  l.O)^  asymptotes 


derived  from  simple  sweep  theory,  may  be  drawn.  The  relative  positions 
of  the  curves  in  figure  5(a)  and  their  asymptotes  suggest  that  a  further 
extension  beyond  x^  will  probably  not  introduce  any  serious  error  in 
these  cases . 

The  curves  are  for  the  most  part  regular  enough  to  permit  interpola¬ 
tion  within  intervals  of  0.2  of  mt-  However,  at  m-^  =  1.0  the  lines 
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diminish  to  a  point  on  the  vertical  axis;  a  curve  for  =0.9  -was 
therefore  inserted  in  the  charts  for  moderate  values  of  m/m^.  With 
m/m-^  less  than  0.5^  m^  =  0.9  represents,  if  x  >xi,  such  extreme 
taper  that  the  successive  points  of  reflection  of  the  Mach  lines 
(xg,  X5, ...)  tahe  place  within  a  very  small  fraction  of  a  chord— length 
and  no  useful  cxirve  can  he  drawn.  No  curves  are  drawn  for  mt  smaller 
than  0.2  because  of  the  tip— interference  limitation  discussed  in  the 
introduction. 

Uncorrected  lift,  region  II.—  In  order  to  find  the  total  lift  (except 
^*0^  tip  losses)  in  region  II,  a  double  integration  with  respect  to  x  and 
y  is  performed  on  equation  (16).  A  first  integration,  with  respect  to 
y,  yields  for  the  indefinite  integral 


qa 


dy  =  4a(x) 


(iiix-Py)(mtCo-mtx+gy) 


+ 


y  mfcCo— (mt-®)x  tan~i 


mx-Py 


(23) 


The  values  of  3y  to  be  sub¬ 
stituted  as  limits  in  equation  (23) 
are  indicated  in  the  sketch  to  the 
right.  Along  the  leading  edge,  the 
right-^iand  member  of  equation  (23) 
reduces  to  zero;  along  the  trailing 
edge  it  becomes 


2jt  cr(x)  Vid^Cq— (m^-^a): 


(24) 


Then  the  total  lift  in  region  II 
(on  both  wing  halves),  except  for  tip 
losses,  is 
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—1  ^2  ^2^3\ 

fi  tan  ^  —  +  — ) 

13  / 


dx  + 


-4. 


CT(x)fidZ  — 


a(x)  tan  ^ 


(25) 


where 


fx 


=  a/  m^c  Q—(  X 


fa  =  V(l+m)(x-xi) 


=  /s/m(x-Ps/m) 


fg  =  y(l+mt)(x2-x)  fg  =  V  lllj.(xij(— x) 


The  indicated  integrations  may  he  performed  n^uIlerically  or  graphi¬ 
cally^  using  values  of  tT(x)  taken  from  the  charts. 

Prinw-ry  tip  correction.—  A  method  of  approximating  the  tip-induced 
reduction  in  loading  is  described  in  reference  4.  The  uncorrected 
loading  along  the  tip  section  is  determined  hy  eqmtion  (I6),  sin^jlified 
hy  replacing  (j(x)  hy  the  constant 


'"s  =  (^)  (26) 

The  assumption  is  then  made  that  the  lift  to  he  canceled  outboard  of  the 
tip  is  a  continuation  of  this  loading  along  lines  parallel  to  the  leading 
edge  of  the  wing.  The  cancellation  is  acconqjlished  hy  means  of  the  '“tip 
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solutions"  of  references  2  and  9> 
flow  fields  having  constant 
pressure  over  a  region  bounded  by 
the  tip  and  a  ray  extending 
outward  from  a  point  Xc^s  on 
the  tip. 

In  canceling  the  cylindri¬ 
cal  field  assumed  in  this  case, 
the  rays  would  all  be  parallel 
to  the  leading  edge.  (See 
sketch. )  The  streamwise  com¬ 
ponent  of  velocity  on  the  wing 
due  to  one  such  element  is 


^c 

-  cos 

at 


m+tQ+2mtj, 

tc^ 


(27) 


where  ug  is  the  constant  value 
of  the  streamwise  perturbation 
velocity  on  the  element  and 
tjj/p  is  the  slope  of  a  ray  from 
its  apex  Xg^s.  This  velocity 
is  multiplied  by  2pV  to 
convert  to  pressure  and  by 


■ 


dS  _  mtg[ct-(^c- 

dtc  p(mt-tc)^ 


to  obtain  the  lift  on  the  differential  of  area  shown  in  the  sketch.  Then 
integrating  with  respect  to  tc  from  —1  to  0  gives  the  total  lift  induced 
on  the  wing  by  the  single  canceling  element  at  xg.  The  result  may  be 
written 


1  dAL  ^ 

- dXp 

qa  dXg 


pVa(  mj.-ffl) 


nUfflS 


mt+m^' 


(1-le)^ 


(28) 


where  Ig,  as  shown  in  the  sketch,  is  the  distance  of  Xg,s  behind  the 
leading  edge,  divided  by  the  tip  chord  c^.,  or 


x^-(ps/m) 


Co+(3s/mt)-Os/m) 


(29) 
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Since  along  the  tip,  from  equation  (I6), 


where  X  is  the  taper  ratio  ct/co>  the  canceling  velocities  Uc  must 
he  given  by  the  differential 


Then  the  primary  tip  coirection  to  the  lift  (both  tips  included) 
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where  K  and  E  are  the  con^jlete  elliptic  integrals  of  the  modulus 


In  the  case  of  an  untapered  wing,  equation  (33a)  reduces  to 


AiL\  — 3jtCo^'/m 

qa  "  2p(l+m)  "'s 


(33b) 


Secondaiy  tip  correction.—  The 
primary  tip  errect  just  aerlved  is 
actually  an  overc6rrection,  since 
the  load  distribution  of  which  it 
is  the  integral  does  not  go  to  zero 
at  the  trailing  edge.  The  residual 
lift  in  the  wake  is  the  result  of 
superposition  of  an  infinite  nuniber 
of  conical  fields  and  therefore 
cannot  be  canceled  identically  in 
any  simple  way.  Numerical  res-ults 
indicate,  however,  as  pointed  out 
in  the  previous  reports  (references 
2  and  4),  that  the  major  part  of 
the  tip  effect  arises  in  the  can¬ 
cellation  of  the  infinite  pressure 
at  the  leading  edge.  In  determin¬ 
ing  the  secondary  tip  correction  to 
the  loading,  therefore,  it  may  be 
assumed  as  an  approximation  that 
the  residual  velocity  field  in  the 
wake  is  conical  with  respect  to  the 
leading--edge  tip,  the  value  of  the 
velocity  along  any  ray  tm  (see 
sketch)  being  detei-mined  by  the  value 


at  the  intersection  x^,y‘^  of  the  ray  with  the  trailing  edge.  Thus  the 
Kutta  condition  is  satisfied,  although  some  pressure  differences  remain 
in  the  wake. 


An  expression  for  the  correction  to  the  total  lift  resulting  from 
this  approximate  cancellation  of  the  tip— induced  velocities  at  the 
trailing  edge  has  been  given  in  reference  1  as  equation  (71),  from  which 
may  be  written 
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(3^) 

where  x*,y*  (see  sketch  p.  15)  is  the  point  of  intersection  of  the  tip 
Mach  line  with  the  trailing  edge,  tm  is  3  times  the  slope  of  a  ray 
from  the  leading-edge  tip,  is  the  intersection  of  the  ray  with 

the  trailing  edge,  Au  is  the  velocity  to  he  canceled  at  x^,yi3, 
and  Au*  is  Au(z*,y*). 

Again  experience  has  shown  the  integral  term  in  eq.uation  (3^)>  which 
involves  very  lengthy  computing,  to  be  considerably  smaller  than  the  term 
in  Au*.  Since  only  an  approximate  tip  correction  is  desired,  it  will 
suffice,  therefore,  to  calculate 

tip 

From  reference  k,  equation  (^5)^ 


Au*  ^  ~^s  (36) 

Ya.  </  mXi* 


where  |*  is  the  streafflwise  distance  of  x*,y*  from  the  leading  edge 
expressed  as  a  fraction  of  the  tip  chordj  that  is. 


i 


(37) 


It  is  easily  determined  that 

3(s-y*)  = 


1+mt 


and  consequently 


mt;(l+m) 

m(  l+nij.) 
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With  these  values,  and  Au*/7a  from  equation  (36),  equation  (35)  becomes 

^  k  gs  i (3A^ 

W  4p  3  1-Hmt  ,/X^  U  “t  ^  ^  ^ 


Total  tip  correction,—  The  total  tip  correction,  to  the  degree  of 
approximation  discussed,  is  the  algebraic  sum  of  equations  (33)  and  (38). 
Except  for  the  occurrence  of  Os^ct  and  \  in  the  coefficient,  the 
resultant  expression  is  a  function  of  m  and  only,  independent  of 


the  tip  location.  Values  of  ^  ^  \  nn  i 

crsct2  /tip 

figure  6,  in  a  form  similar  to  the  chart  of 


have  been  plotted  in 

(fe)^ 


Numerical  Example 


As  a  summary  of  the  method,  a  sample  calculation  will  be  outlined. 
The  lift— curve  slope  will  be  calculated  for  an  untapered  wing  of 

10— foot  chord  and  40-foot  span,  swept  back  45®  and  flying  at  a  Mach 
number  of  I.08.  Then  cot  A  =  1.0,  p  =  0.4,  m  =  0.4,  mt  =  0.4  and, 
from  figure  4, 


L  ^  _ _ m^ 


qa 


=  13*30--^  Cq^  =  1330  square  feet 


I  3 


(it  should  be  noted  that  ^  =  cot  a) 

P 

With  values  of  a  obtained  from  figure  5(a) ^  graphical  integration  of 
equation  (25)  gives 


f-) 

Vqa 


=  531  square  feet 


Since  the  leading--edge  sweep  is  45®,  the  tip  is  at  x  =  20,  and,  from 


figure  5(a),  ag=0.635(  “  O.518.  Erom  figure  6,  since  A  “  1, 


2 

AT  N 

—  )  =  — 1.69  — r —  =  —219  square  feet 


tip 


3 
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Combining  *the  three  components  of  the  lift^  we  obtain  -  1642  sciuare 
feet  and 

ERAG  DUE  TO  LIFT 


The  drag  due  to  lift  can  be  found  as  the  integral  around  the  airfoil 
of  the  incremental  pressure  puV  times  the  slope  w/v*  of  the  surface. 

In  the  case  of  a  thin  flat  airfoil,  this  calculation  gives  merely  the 
lift  times  the  angle  of  attack  over  most  of  the  surface.  However,  as  is 
discussed  in  reference  11  in  connection  with  the  two-dimensional  wing  in 
subsonic  flow,  an  infinite  pressure  acting  on  the  leading  edge  results  in 
a  suction  force,  which  tends  to  reduce  the  total  drag.  When  the  wing  is 
swept  behind  the  Mach  lines  from  its  apex,  a  similar  suction  can  exist  in 
supersonic  flow,  since  the  pressure  on  a  subsonic  leading  edge  is  theoret¬ 
ically  Infinite. 

Hayes  (references  7  12),  Robinson  (see  reference  8),  and  others 

have  derived  the  formula  for  the  suction  force  on  a  subsonic  leading 
edge  by  assuming  the  flow  near  the  leading  edge  to  be  essentially  two- 
dimensional  and  applying  the  results  of  two-dimensional  potential  theory. 
The  single  result  obtained  in  that  manner  has  been  verified  for  the  long 
swept— back  wing  by  application  of  the  somewhat  different  approach  of 
reference  I3. 

By  the  two-dimensional  approach,  the  suction  force  is  foimd  to  be 
proportional  to  the  square  of  the  strength  of  the  leading-edge  singularity 
in  the  perturbation  velocity  u.  This  quantity  has  already  been  discussed 
in  connection  with  the  lift  in  region  II,  where  the  strength  of  the 
leading-edge  singularity  was  defined  as  the  coefficient  of  the  inverse 
sqiiare  root  of  mx— Py.  This  coefficient  is  (equation  (19))  the 

forward  part  of  the  wing  and  Ca  +  (^)o  ^ 

o 

behind  the  -  tralling-edge  Mach  line.  From  Robinson’s  work,  the  longitudi¬ 
nal  component  of  the  suction  force  per  unit  length  in  the  x  direction 
may  be  expressed,  for  region  I,  in  terms  of  as  follows: 


dx 


pn 


m 


and  similarly  for  the  remainder  of  the  wing,  with  replaced  by  the 

corrected  coefficient. 
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For  a  given  wing,  the  total  drag  due  to  lift  is  the  result  of  sub¬ 
tracting  the  total  thrust  from  the  product  of  the  lift  and  the  angle  of 
attack;  that  is> 


D 


dT 

dx 


dx 


or,  in  coefficient  form. 


(39) 


(40) 


where  S  is  the  total  area  of  the  wing.  Writing 


we  obtain 


2 


3s 

dx 


=  It 


(42) 


so  that 


(43) 


with  the  final  integration  to  be  performed  numerically  or  graphically. 
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APELICATIOK  AM)  DISCUSSION 


Lift-Curve  Slope 


The  lift-curve  slope  has  heen  calculated  for  two  families 

of  untapered  wings ^  with  m  constant  at  0.2  and  O.k,  and  with  varying 
aspect  ratio.  The  results  are  plotted  against  the  reduced  aspect  ratio 
pA  in  figure  7*  The  circled  points  indicate  the  actual  cases  for  which 
calcinations  were  made;  the  faired  curves  are  only  approximations,  since 
they  ignore  slight  discontinuities  in  curvature  associated  with  the  onset 
of  interaction  between  the  wake  and  leading  edge,  and  higher— order  dis¬ 
continuities  associated  with  successive  reflections  of  the  trailing— edge 
Mach  lines. 

Two  points  corresponding  to  tapered  wings  with  m  =  0.4  and 
iQ^  =0.6  are  also  Included.  In  one  case  (pA  =  1*6),  and  in  the  case 
of  the  untapered  wing  with  the  same  span  (PA  =  1.2),  the  trailing-edge 
Mach  lines  did  not  intersect  the  leading— edge ,  and  the  values  of  Cj^ 
were  obtained  by  the  formulas  of  reference  1.  It  should  also  be  men¬ 
tioned  that  the  remaining  values  agree  within  2  or  3  percent  with  values 
calculated  entirely  by  the  conical-flows  method,  the  forn^r  being  slightly 
lower  than  the  latter. 

The  curves  for  the  untanered  wings  may  be  seen  to  be  approaching  at 
the  upper  end  the  value  23(m/  ^1— given  by  simple  sweep  theory.  At 
the  lower  end,  the  curves  should  approach  the  origin  along  the  line 

Ct  =  -  a  given  by  low-aspect— ratio  theory  (reference  l4).  As  previously 
■MX  2 

mentioned,  the  present  calculations  cannot  properly  be  extended  below 
pA  =  1  because  of  Interference  between  the  flow  fields  originating  at  the 
tips.  However,  two  such  cases  have  been  included  for  m  =  0.2  because, 
with  so  much  sweep,  the  wing  areas  affected  are  small  and  the  interference 
losses  should  be  negligible.  Results  in  these  latter  cases  may  be  compared, 
because  of  the  small  apex  angle  relative  to  the  Mach  angle,  with  the  results 
of  the  slender^ing  theory  of  reference  5.  This  comparison  is  shown  in 
figure  7,  although  a  discrepancy  in  plan  form  lessens  the  significance  of 
the  agreement.®  The  results  of  reference  5  have  also  been  plotted  for 
m  =  0.4,  in  which  case  the  assun^tion  of  extreme  slenderness  is  no  longer 
justified  and  introduces  an  appreciable  error.  (It  should  be  mentioned 
that  the  asymptote  for  the  slender-wing-theoryjcurves  is  below  the  value 
given  by  simple  sweep  theory  by  the  factor  vl-m^,) 

®The  theory  of  reference  5>  while  applicable  to  any  swept— back  wing  lying 
well  within  the  leading— edge  Mach  cone,  has  been  worked  out  only  for  a 
liiaited  family  of  plan  forms,  having  straight  leading  and  filleted 
trailing  edges,  and,  consequently,  a  slight  taper. 
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Drag  Due  to  Lift 


In  figure  8,  the  drag-rise  factor  Cd/Cl^  is  plotted  for  the  same 
families  of  vings.  Comparison  is  made  with  a  theoretical  minimum  for 
slender  wings  in  supersonic  flight  obtained  by  R.  T.  Jones  in  an  unpub— 
lished  analysis.  Using  a  method  similar  to  Hayes*  (reference  12)  and 
assuming  the  wing  to  be  narrow  compared  with  the  Mach  cone ,  Jones  has 
derived  a  minimum  "wave  drag"  coefficient 


2rtA^ 


where  is  the  aspect  ratio  defined  in  the  streamwise,  instead  of  the 

spanwise,  direction;  that  is^  if  Z  (numerically  equal  to  x^)  is  the 
ove]>-all  length  of  the  wing. 


m 


The  wave  drag  is  to  be  added  to  the  "vortex  drag,"  which  is  the 
induced  drag  of  subsonic  flow,  calculated  from  the  spanwise  loading. 
Using  the  minimum  induced  drag  calculated  from  lifting-line  theory  gives 
as  the  minimum  supersonic  drag-rise  factor 

-Si  =  X  +  _&f_  (46) 

Cjf  itA  2ttAx 


It  may  be  seen  that  the  drag  rise  of  the  constant— chord  swept-back 
wings  is  fairly  close  to  the  minimum,  especially  at  the  lower  values  of 
m  for  which  equation  (46)  was  derived. 


Ames  Aeronautical  Laboratory, 

National  Advisory  Committee  for  Aeronautics, 
Moffett  Field,  Calif.,  Mar.  15,  1950. 
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APEENDIX 

SYSTEM  OF  NOTATION 


General 


free-stream  velocity 
free-etream  Mach  number 

Ju^-1 

density  of  air 

/ 1  2^ 

dynamic  pres  star  e  I  -^pV  ) 

pressiore  difference  between  upper  and  lower  surfaces^  or 
local  lift 

angle  of  attack,  radians 


lift 

lift  coefficient 
drag  coefficient 


Wing  Dimensions 


root  chord 
tip  chord 
semispan 
wing  area 

over— all  length  of  the  wing  in  the  streamwise  direction 
angle  of  sweep  of  the  leading  edge 
taper  ratio  (c^/cq) 
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A  aspect  ratio  (4s2/s) 

Coordinates 

Xyj  Cartesian  coordinates  in  the  stream  direction  and  across  the 
stream  in  the  plane  of  the  wing 

^^^a  coordinates  of  apex  of  oblique  trailing-edge  eleTOnt  (see 
equation  (l4)) 

^b^y^b  coordinates  of  point  on  trailing  edge^  within  the  tip  Mach  cone 

x^^s  coordinates  of  point  on  tip;  apex  of  tip  element 

coordinates  of  intersection  of  trailing-edge  Mach  cone  with 
leading  edge  (see  equation  (ll)) 

coordinates  of  intersection  of  Mach  line  from  with 

trailing  edge  (see  equation  (13)) 

coordinates  of  intersection  of  tip  Mach  line  with  trailing  edge 

Xrp^s  coordinates  of  intersection  of  tip  and  trailing  edge 

I*  streamwise  distance  of  x*^y^  back  from  leading  edge,  as  a 

fraction  of  the  tip  chord  (equation  (37)) 

distance  of  x^  behind  leading-edge  tip,  as  a  fraction  of  the 
tip  chord  (equation  (29)) 


In  the  following,  all  slopes  are  measured  counterclockwise  from  a 
line  extending  downstream  from  the  apex  of  the  wing  or  of  the  pertinent 
elementary  sector: 


m 


“t 

a 


slope  of  leading  edge  _  ^ 
slope  of  Mhch  lines 

slope  of  trailing  edge 
slope  of  Mach  lines 


slope  of  ray  from  the  origin 
slope  of  Mach  lines 


slope  of  ray  from  trailing-edge  apex 
slope  of  Mach  lines 


X-Co 
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slope  of  ray  from  leading-edge  tip  y-s 

■  I  -.1  .1.1  ■  I  ■  I  .  . .  .1. 1..—— . .  r:  p  ■  ■' 

slope  of  Mach  lines  x—  ^ 

m 


slope  of  ray  from  apex  of  element  a 
slope  of  Mach  lines 

slope  of  ray  from  apex  of  tip  element 
slope  of  Mach  lines 


y-ya 

X-Xq 


the  value  of  a  corresponding  to  a  trailing-edge  element  of 
which  the  apex  lies  on  the  Mach  forecone  of  the  point  at  which 
the  loeui  is  being  calculated  (equation  (22)) 

(equation  (8)) 

the  smelliest  value  of  to  along  which  cancellation  of  pressure 
ahead  of  the  leading  edge  can  affect  the  point  at  which  pressin-e 
is  being  calculated.  For  a  point  on  the  leading  edge, 

To  =  "to  =  iiix/(2:-Co) 

the  smallest  values  of  tg^  along  which  cancellation  of  pressure 
ahead  of  the  leading  edge  can  affect  the  point  at  which  pressure 
is  being  calculated.  For  a  point  on  the  leading  edge, 

^a  =  ^a  =  (inx-(3yg^)/(x-Xg^) 


Streamwise  CoiJ5)onents  of  Perturbation  Velocity 

basic  (uncorrected)  ijertiirbation  velocity  as  given  by  solution 
for  triangular  wing  (equation  (l)) 

value  of  Uj^  at  a=0  (equation  (2)) 

constant  perturbation  velocity  on  canceling  (oblique)  sector  in 
wake  (equation  (6)) 

constant  perturbation  velocity  on  canceling  sector  outboard 
of  tip 

symmetrical  trailing-edge  correction  to  u^  (equation  (3)) 

da,  correction  to  ua  due  to  single  oblique  trailing-edge 
da 

element  (equation  (5)) 

value  of  tip  correction  to  u  at  the  point  x*^y*  (equation  (36)) 
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Artitrary  Mathematical  Symbols 

^“1  /  2 

value  of  coefficient  of  (iiixH3y)  in  at  the  leading 

edge  ( equat ion  ( 19 ) ) 

(AC)q  decrement  in  due  to  reflection  of  (Au)q  at  leading  edge 

(equation  (20;) 

da  decrement  in  due  to  reflection  of  (Au)q^  at  leading  edge 

(equation  (21)) 

o  leading-edge  correction  coefficient  defined  hy  equation  (I7) 

Oq  value  of  a  at  leading-edge  tip 

Elliptic  Integrals 

modulus  of  elliptic  integral,  defined  where  used  (also  with 
subscripts) 

argument  of  elliptic  integrals,  defined  where  used  (also  with 
subscripts) 

incomplete  elliptic  integral  of  the  first  kind  of  modulus  k 
and  argument  cp 

complete  elliptic  integral  of  the  first  kind  of  modulus  k; 
that  is,  K  =  F(k, 

iiicoii5)lete  elliptic  integral  of  the  second  kind 

contplete  elliptic  integral  of  the  second  kind  of  modulus  k; 
E(k,|) 

BEFEREIirCES 


1.  Lagerstrom,  P.  A.,  Wal3.,  D.  D.,  and  Graham,  M.  E.:  Formulas  in  Three- 

Dimensional  Wing  Theory.  Douglas  Eep.  SM-II9OI,  Aug.  19^6. 

2.  Cohen,  Doris:  The  Theoretical  Lift  of  Flat  Swept— Back  Wings  at 

Supersonic  Speeds.  HACA  TN  1555#  19^8. 


k 

cp  or  >|f 

F(k,«p) 

K,K(k) 

E(k,(p) 

E,E(k) 


26 


NACA  TN  2093 


3.  Malvestuto,  Frank  S.,  Jr.,  Margolis,  Kenneth,  and  Eihner,  H.  S.: 

Theoretical  Lift  emd  Damping  in  Eoll  of  Thin  Swepthack  Wings  of 
Arbitrary  Taper  and  Sweep  at.  Supersonic  Speeds.  Subsonic 
Leading  Edges  and  Supersonic  Trailing  Edges.  HAGA  TN  i860,  19^9- 

4.  Cohen,  Doris:  Theoretical  Loading  at  Supersonic  Speeds  of  Flat 

Swept-Back  Wings  with  Interacting  Trailing  and  Leading  Edges . 

NACA  TN  1991,  1949. 

5.  Lomax,  Harvard,  and  Heaslet,  Max.  A.:  Linearized  Lifting-Surface 

Theory  for  Swept-HBack  Wings  with  Slender  Plan  Forms.  NACA 
TN  1992,  19^9. 

6.  Jones,  R.  T.:  Estimated  Lift-Drag  Ratios  at  Supersonic  Speed. 

NACA  TN  1350,  19^7. 

7.  Hayes,  W.  D.,  Browne,  S.  H.,  and  Lew,  E.  J.:  Linearized  Theory  of 

Conical  Supersonic  Flow  with  Application  to  Triangular  Wings, 

North  American  Aviation  Rep.  NA— k6— 818,  Sept.  30,  19^6. 

8.  Squire,  E.  B.:  An  Example  in  Wing  Theory  at  Supersonic  Speeds. 

Abstract  in  Proceedings  of  the  Seventh  International  Congress 
for  Applied  Mechanics,  19^8.  vol.  2,  Part  II,  p.  6l4. 

9.  Lagerstrom,  P,  A.:  Linearized  Supersonic  Theory  of  Conical  Wings. 

NACA  TN  1685,  1948. 

10.  Busemann,  A.:  Infinitesimal  Conical  Supersonic  Flow.  NACA 

TM  1100,  1947- 

11.  Betz,  A.:  Applied  Airfoil  Iheory.  Vol,  IV,  dlv.  J  of  Aerodynamic 

Theory,  W.  F.  Durand,  ed.,  Julius  Springer  (Berlin)  1935*  P*  28. 

12.  Hayes,  W.  D.:  Linearized  Supersonic  Flow.  North  American  Aviation, 

Inc.,  Rep.  no.  AL-222,  June  18,  1947* 

13.  Jones,  E.  T.:  Leading-Edge  Singularities  in  Thin  Airfoil  Theory. 

JouTo  Aero.  Sci.,  vol.  17,  no.  5,  May  1950,  pp.  307-310. 

14.  Jones,  E.  T.:  Properties  of  Low-A.spect-Eatio  Pointed  Wings  at  Speeds 

Below  and  Above  the  Speed  of  Sound.  NACA  Rep.  835,  1946. 


Mach  number 


NACA  TN  2093 


Figure  /.  -  Plan  view  of  swept-bock  wing  at  two  supersonic 
Mach  numbers. 


Figure  2.-  Range  of  applicability  af  farmutas  and  charts. 
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Figure  4.-  Chart  for  the  computation  of  the  lift  in  region  X . 
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Figure  5.-  Charts  for  determining  <r,  the  strength  of  the  leading-edge  singularity. 
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Figure  6.-  Chart  for  the  correction  of  the  lift  for  tip  effect. 


